Topological spin Hall and spin Nernst effects in a bilayer graphene 
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We consider intrinsic contributions to the spin Hall and spin Nernst effects in a bilayer graphene. 
The relevant electronic spectrum is obtained from the tight binding Hamiltonian, which also includes 
the intrinsic spin-orbit interaction. The corresponding spin Hall and spin Nernst conductivities are 
compared with those obtained from effective Hamiltonians appropriate for states in the vicinity of 
the Fermi level of a neutral bilayer graphene. Both conductivities are determined within the linear 
response theory and Green function formalism. The influence of an external voltage between the 
two atomic sheets is also included. We found transition from the topological spin Hall insulator 
phase at low voltages to conventional insulator phase at larger voltages. 

PACS numbers: 73.43.-f, 72.25. Hg, 73.61.Wp 



I. INTRODUCTION 

Four decades ago Dyakonov and Perel showed that 
a system with strong spin-orbit interaction should re- 
veal transverse spin current and spin accumulation in the 
presence of external longitudinal electric fieldi 2 - - even 
if the system is nonmagnetic. This effect, known now as 
the spin Hall effect (SHE)£, was studied extensively in the 
last few years^r— , and is still of current interest - mainly 
because it offers a new possibility of spin manipulation 
with electric field only. The possibility of pure electrical 
manipulation of spin degrees of freedom is interesting not 
only from fundamental reasons, but also from the point of 
view of possible applications in future spintronics devices 
and information processing technologies^—. 

The crucial interaction responsible for SHE, i.e. the 
spin-orbit coupling, may be either of intrinsic (internal) 
or extrinsic origin. The corresponding extrinsic SHE is 
associated with mechanisms of spin-orbit scattering on 
impurities and other defects (skew scattering and/or side 
jump), while the intrinsic SHE is a consequence of a 
nontrivial trajectory of charge carriers in the momen- 
tum space due to the spin-orbit contribution of a perfect 
crystal lattice to the corresponding band structure. The 
intrinsic SHE may be described in terms of the Berry 
phase formalis m 14 : 15 and therefore it is also referred to 
as the topological SHE. 

It is well known that various spin effects, like for in- 
stance spin current and spin accumulation, may be gen- 
erated not only by external electric field, but also due to 
a temperature gradient. Indeed, there is a great interest 
currently in spin related thermoelectric effects. One of 
such phenomena is the spin Sccbcck effect, where longi- 
tudinal spin current and spin voltage are generated by a 
temperature gradient^. Of particular interest, however, 
are spin thermoelectric effects in systems with spin-orbit 
interaction, where a temperature gradient gives rise to 
transverse spin accumulation and/or spin currents. Thus, 
the temperature gradient in such systems may lead to 
anomalous (in case of ferromagnetic systems) and spin 



FIG. 1. (color on-line) Schematic presentation of transverse 
spin accumulation induced by (a) electric field (spin Hall ef- 
fect) and (b) temperature gradient (spin Nernst effect). 



Nernst effects^—, which correspond to the anomalous 
and spin Hall effects induced by longitudinal electric field. 
Similarity of the spin Hall and spin Nernst effects in non- 
magnetic systems with spin-orbit interaction is presented 
in Fig.l, which clearly shows that the SHE is generated 
by external electric field, while the spin Nernst effect 
(SNE) is a similar effect generated by a temperature gra- 
dient instead of electric field (gradient of electrostatic 
potential). 

In this paper we consider the topological contribu- 
tion to the spin Hall and spin Nernst effects in a bi- 
layer graphene. Graphene is a two-dimensional crystal of 
carbon atoms. A monolayer graphene in a free standing 
form was obtained few years ago and owing to its unusual 
and peculiar properties quickly became one of the most 
extensively studied materials 2 ^—. After the pioneering 
paper by Kane and Mele 2 ^, spin Hall effect in a mono- 
layer graphene was studied in many papers and in various 
physical situations. Recently, bilayer graphene is exten- 
sively studied as a more appropriate for applications than 
a single-layer one. Moreover, it has been shown that spin- 
orbit coupling in a bilayer can be enhanced in comparison 
to that in a single-layer graphene 2 ^. This motivated us to 
consider topological contributions to the SHE and SNE 
in a bilayer graphene. 

Both SHE and SNE in (nonmagnetic) graphene are 
generated by spin-orbit interaction. In general, one can 
distinguish two different forms of the spin-orbit interac- 
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tion having the crystal lattice periodicity and contribut- 
ing to the relevant band structure - internal and Rashba 
spin-orbit interactions. The latter interaction is due to a 
substrate and can be controlled by an external gate volt- 
age. In this paper we consider the contributions to SHE 
and SNE due to the intrinsic spin-orbit interaction only. 
It is known that this interaction opens an energy gap 
at the Dirac points. It was also shown that the energy 
gap can be tuned externally by applying a voltage bias 
between the layer o 25 ' 26 . The gate voltage dependence of 
the spin Hall conductivity, leading to phase transition be- 
tween the spin Hall insulator and classical insulator will 
also be considered in this paper. 

The description of transport properties of graphene is 
usually based on some effective Hamiltonian, which prop- 
erly describes electronic spectrum near the Fermi level of 
a neutral system. However, it is well known that the 
topological contribution to the spin Hall effect includes 
contributions from electronic states far below the Fermi 
level, and therefore a more accurate electronic spectrum 
is required to describe the effect properly. Accordingly, 
in this paper we calculate the spin Hall and spin Nernst 
conductivities from a more realistic electronic spectrum 
based on a tight binding Hamiltonian, and compare them 
with the corresponding conductivities obtained on the 
basis of effective Hamiltonians. The results presented in 
this paper reveal, however, a very good agreement be- 
tween the conductivities derived from the effective and 
tight binding Hamiltonians. 

The paper is organized as follows. In section 2 we de- 
scribe briefly the electronic states of a bilayer graphene 
within the tight binding Hamiltonian and also in terms 
of an effective Hamiltonian, which is sufficient for states 
in the vicinity of the Fermi level of a neutral graphene. 
In both cases the intrinsic spin-orbit interaction and the 
effect of a normal bias voltage are taken into account. In 
section 3 we calculate the spin Hall and spin Nernst con- 
ductivities for the tight binding Hamiltonian and com- 
pare them with those obtained from the effective Hamil- 
tonian. In the latter case we derive analytical formulas 
for the spin Hall and spin Nernst conductivities. We also 
discuss the role of normal bias (in the framework of the 
effective model). Description based on a reduced low- 
energy effective Hamiltonian is presented in section 4. 
Summary and final conclusions arc given in section 5. 



II. ELECTRONIC SPECTRUM OF THE 
BILAYER GRAPHENE 

A single-layer graphene is a monolayer of carbon 
atoms arranged in a two-dimensional honeycomb lat- 
tice which can be also considered as being composed of 
two nonequivalent triangular sublattices. In the absence 
of spin-orbit interaction the Fermi surface of a neutral 
single-layer graphene consists of two nonequivalent K 
and K' points of the Brillouin zone, at which the va- 
lence and conduction bands touch each other. The cor- 



responding electronic spectrum can be described by a 
tight binding Hamiltonian with nearest and next-nearest 
neighbor hopping terms. The low energy electron states 
near the points K and K' can be well approximated by a 
conical energy spectrum (linear dispersion relations) . As 
a result, charge carriers in the vicinity of the points K 
and K' arc described effectively by the relativistic Dirac 



equation 
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Intrinsic spin-orbit interaction opens then 



an energy gap at the Dirac points^. The tight bind- 
ing and effective Hamiltonians for a bilayer graphene are 
more complex, as described below. 



A. Tight binding model 

The bilayer graphene in the Bernal stacking (see eg 
[22T |) is described by the following tight binding Hamilto- 
nian: 



where 



H 



H = 



H r 



d 2 kV(k)T I rt H J V(k) 



h so S z + VSq hoSo 



(1) 



h*S -h so S z -VS )' (2) 



r = 



71 S 




(3) 



The matrix elements ho and h so are defined as fol- 
lows: ho = -t[e abk y/' i + 2cos{\ak x )e- lbk y/ 3 ] and h so = 
— 2i'[sin (ak x ) — 2s'm(ak x /2) cos (bk y )], where t is the 
hopping integral between the nearest neighbors in the 
atomic sheets, t' is the next-nearest neighbor spin-orbit 
hopping amplitude, while b = aV3/2 with a being the lat- 
tice parameter. Furthermore, V is the voltage between 
the two atomic sheets of the bilayer (measured in energy 
units), S a denote the unit (a = 0) and Pauli (a — x, y, z) 
matrices in the spin space, while 71 describes coupling 
between the two atomic layers. 

The corresponding energy eigenvalues for V = have 
then the following form: 



or 



2hl 



2h% 



Ti 



1 1/2 



(4) 



-1 1/2 



(5) 



This spectrum is shown in Fig. 2(a). States near the point 
K are shown by the solid lines in parts (b) and (c). The 
part (c) reveals a small energy gap created at the Dirac 
points by the spin-orbit interaction. 

When V 7^ 0, the inversion symmetry is broken (layers 
are no longer equivalent) and the degeneracy is lifted. 
The corresponding eigenvalues acquire then the form 



Ei. 



hi + hi + V 2 + ?l 



i[(7 1 2 -4yM 2 + 4( 7l 2 + 4^)/ l 2 ] 1 / 2 



1/2 



,(6) 
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FIG. 2. Electronic spectrum of the bilayer graphene in the 
tight binding model for k v = 0, t = 2.9 eV, t' /t = 1CT 5 , 
7i/t = 0.069, and V — (a). The tight binding energy 
spectrum in the vicinity of the K point (solid lines) and the 
corresponding energy spectrum obtained from the effective 
model (dotted lines) are compared in (b) for V = and in 
(d) for V/t — 0.01. Parameters of the effective model, cor- 
responding to the tight binding one, are: A so = 0.15 meV 
and v — 3.516 x 10~ 10 eVm. The electron states very close 
to the energy gap are shown in (c) for V — and in (c) 
for V/t — 0.01 (here the results from the tight binding and 
effective Hamiltonians overlap). 



Ei> 2 > = T 



2 . 7? 



1/2 



(7) 



E3A — T 



hi + h* + V 



2 , 7i 
2 



1/2 



,(8) 



E31 A> — T 



K + h 2 so + V 2 + 4 



i[( 7l 2 + 4 yM 2 + 4(7 1 2 + 4y 2 )^] 1 / 2 



1/2 



(9) 



The corresponding spectrum for the assumed value of 
V/t = 0.01 is indistinguishable from the spectrum shown 
in Fig. 2(a) for V = 0. The differences can be seen on a 
smaller energy scale, as in the parts (d) and (e) of Fig. 2. 
When comparing Figs 2(b) and 2(d), one can notice a 
larger energy gap for V/t = 0.01. This is more clearly 
seen when comparing the parts (c) and (c). First, for 
the assumed value of V the gap is wider than for V = 0. 
Second, the top and bottom band edges become split and 
shifted away from the Dirac points. As will be described 
later, the applied voltage between the atomic sheets first 
closes the gap and then opens a new one with the width 
increasing with V. 



B. Effective Hamiltonian 

When only electronic states near the Fermi level (near 
the Dirac points) of a neutral graphene are relevant, one 
can make use of some effective Hamiltonians to describe 
the corresponding electronic spectrum. Such a Hamilto- 
nian can be derived using the k • p approximation. As 
a result, the effective Hamiltonian for states near the K 
point of the bilayer graphene takes the for m 27 ' 28 : 



h k = %® m 



.21(7; 
2 [ ' x 



o x ® S - T y ® <Jy ® S ), 



(10) 



where 



H S K = v(k x a x + k y (Ty) ® #0 

+A so a z ®S Z + Voo ® So- (11) 

The first term on the right side corresponds to two de- 
coupled atomic monolayers, each of them being described 
by the Kane Hamiltonian for a single layer graphene, 
H S K . In turn, the second term describes coupling be- 
tween the monolayers, with T a denoting respectively the 
unit matrix (a = 0) and Pauli (a = x,y,z) matrices 
associated with the layer degree of freedom. In turn, 
a a (a = 0, x, y, z) are the unit (a = 0) and Pauli 
(a = x, y, z) matrices in the pseudo-spin (sublatticc) 
space. Relations between the parameters of the tight 
binding and Kane models are: v = sfita/l = hvp {vf is 
the carrier velocity at the Fermi level) and A so = 3%/3t'. 

The eigenvalues of Hamiltonian ([T7J| for V = take 
the form 



E h2 = T 
and 

E3A = T 



fcV + ^+A 2 -^v/4fc 2 z; 2 + 7 2 



,1/2 



(12) 



fcV + 21 + Al + 21^/W + 7i 2 



1/2 



(13) 
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Electronic spectrum near the point K, described by the 
above formula, is shown in Fig.2(b,c), where it is com- 
pared with the spectrum obtained from the full tight 
binding Hamiltonian. Close to the K point (gap), spectra 
from both models coincide very well. 
When V ^ 0, one finds 



£1,2 = T 



v 2 k 2 + Ai 



if 



i[( 7l 2 - WAJ 2 + 4 v 2 fc 2 ( 7l 2 + 4V 2 )] 1 / 2 



1/2 



where cr?j > (for i,j = x,y) is the spin Hall conductiv- 
ity with s n — (h/2)a n being the 77-th component (77 = 
x, y, z) of electron spin, while a?" denotes the thermoelec- 
tric spin Nernst conductivity. The two conductivities are 
not independent and obey some general relations. Our 
objective is to hnd first the zero-temperature spin Hall 
conductivity, and then to calculate the low-temperature 
thermoelectric spin Nernst conductivity from these rela- 
tions, as described below. 

The quantum-mechanical operator of spin current den- 
sity may be defined as 



Ev,2' = T 
f 



'Y 2 

v 2 k 2 + A 2 + V 2 + 



[( 7l 2 + 4VA so ) 2 + 4« 2 fc 2 ( 7l 2 + 4V 2 )] 1 / 2 



-1 1/2 



05) 



£3,4 = T 



v 2 k 2 + A 2 n + V 2 + i 



i[( 7l 2 - WA so ) 2 + 4<; 2 fc 2 ( 7 2 + 4V 2 )] 1 / 2 



nl/2 



E, 



v 2 k 2 + A 2 n 



V 



2 , 7i 2 



+ i[( 7l 2 + WA so ) 2 + 47; 2 fc 2 ( 7l 2 + 4V 2 )] 1 / 2 



(16) 



1/2 
(17) 



The above spectrum is shown in Fig.3(d,e), where it is 
compared with the corresponding spectrum obtained in 
the tight binding model. As before, spectra from tight- 
binding and effective models coincide near the K point. 
Note, the band splitting due to V is well resolved only in 
part (e). 

Separation of the bands E% and E4 in the effective 
model described above, as well as in the tight binding 
model, is much larger than the separation of the bands 
Ei and E2, see Fig. 2. Therefore, when the electronic 
states close to the band edges are relevant and sufficient 
to describe transport properties (eg. when the Fermi 
level is in the gap), one may restrict considerations to the 
bands E\ and E^. This leads to a further simplification 
of the effective Hamiltonian, as described in more details 
in section 4. 



III. SPIN HALL AND SPIN NERNST EFFECTS 

Spin Hall and spin Nernst effects correspond to 
transversal spin currents induced by electric field and 
temperature gradient, respectively. By analogy to the 
usual Hall and Nernst effects one may write the density 
of spin current due to electric field E and temperature 
gradient VT as 



(18) 



[V,S„]_ 



(19) 



where [A, B] + = AB + BA denotes the anticommutator 
of any two operators A and B, while u, = (l/K)(dH/dki) 
is the velocity operator (i = x,y). The latter operator 
can be easily found from the corresponding Hamiltonian 
[Eqs (1) and (10)]. In the linear response theory, the 
frequency-dependent spin Hall conductivity is then given 
by the formula^, 

= £ Tr / £(2^ K^]+Gk(e + c,KG k ( £ ), 

(20) 

where G^e) is the Green function corresponding to the 
appropriate Hamiltonian of the system. When we re- 
strict considerations to the topological contribution to 
the spin Hall current in the d.c. limit, this formula gives 
exactly the same result as that based on the Berry phase 
calculations^ in momentum space. 

It has been shown that the Berry phase leads to an 
additional term in the general expression for the orbital 
magnetization^. This correction gives rise to some con- 
tributions to the charge and spin current o 17 ' 31 , and also 
allows to write the relationship between intrinsic spin 
Nernst conductivity and intrinsic zero-temperature spin 
Hall conductivity, which in the low-temperature regime 
takes the form 3 -^ 



3e 



T 



ds 



(21) 



where T stands for temperature and ks denotes the 
Boltzman constant. The latter equation is the spin ana- 
log of the Mott relation for charge transport, and will be 
used to calculate the low-temperature spin Nernst con- 
ductivity from the zero-temperature spin Hall conductiv- 
ity. The derivative in Eq.(21) is taken at the Fermi level 
li. The latter can be tuned by an external gate voltage. 

Thus, we need to calculate the spin Hall conductivity 
first. The relevant derivation depends on the model ap- 
plied to describe the corresponding electronic spectrum. 
Below we present derivation of the conductivity for the 
effective Hamiltonian, where analytical results are avail- 
able. These results will be compared with those obtained 
numerically for the tight binding model. 
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A. The limit of V = 

Assume first the limit of V = 0. To find the spin Hall 
conductivity we start from Eq. (|20[) and write it in the 
form 

e f de f d 2 k 



xyK ' 2co J 2n J (2tt) 2 ' 

4 



X 



~[[e- E n +us + fi + iS sign(e)] 2 

4 

x JJ [e - £ m + M + iS Bign(e)] ~ 2 . (22) 

m— 1 

Here, D(e + w,e) is defined as 

L>(e + w,e) =Tr{[u a!) s a ] + 5k(£: + a;)u y 5k(e)}, (23) 

where gu_(e) denotes the nominator of the corresponding 
Green function Gk(e). Taking the first two terms of the 
expansion of D[e + uj, e) with respect to oj, one finds 

D( E + u,e)2iiwx(e), (24) 

with 

X (e) = 8v 2 A so {v 2 k 2 [v 2 k 2 + 2(A 2 - (s + M ) 2 )] 

+(A 2 -( £ + ^) 2 )(7i 2 +A 2 -(£ + ^) 2 )} 2 
x{(A 2 -( e + ^) 2 )( 7l 2 +A 2 -( £ + / ,) 2 ) 
+ , 2 fc 2 [ u 2 fc 2 + 2( 7l 2 +A 2 -(e + M) 2 )]}.(25) 

Thus, in the limit of uj — > one finds the following ex- 
pression for the spin Hall conductivity 



J 2nJ (2^^ h 



where 



Integrating now over e one finds 



(26) 



(27) 



(28) 



where R n (n = 1 — 4) are the residua associated with the 
corresponding selfenergies (electron bands), and f(E) is 
the Fermi distribution function (here for zero tempera- 
ture). These residua are equal: 



Rl,2 = ± 

and 



4V2v 2 A so [2v 2 k 2 ( 7l + + 7i (7? - 7i£ + 2A 2 J] 



^3,4 — ±- 



£3 (2 w 2 fc 2 + 7l 2 - 7l £ + 2A 2 ) 3/2 
8y/2v 2 A so L 



£3 (2v 2 k 2 + 7l 2 + 7 i£ + 2A 2 Q ) 



5/2 



(29) 
(30) 




-0.4 -0.2 0.2 0.4 
H [meV] 

FIG. 3. (color on-line) Spin Hall and spin Nernst conductiv- 
ities of a bilayer graphene for V = and other parameters 
as in Fig. 2. Contributions from both K and K 1 points are 
included. The solid lines correspond to the results obtained 
from the effective Hamiltonian (10), while the dotted ones 
from the tight binding model. 



with 



L = - v 2 fc 2 (3 7l 3 + 7l 2 £ + 4 7l A 2 Q - 2£A 2 J 

+2w 4 fc 4 ( 7l +0 
-7i (7i 4 + 7i£ + 2 7l 2 A 2 Q + 2 7l £A 2 Q + 2A 4 D ), (31) 



where £ = \J Av 2 k 2 + 7 2 . Thus, the conductivity may be 
written in the form 



°% = -^E / kR n f(E n )dk. 



(32) 



Taking into account the integrals 

£(2, 2 P+ 7l 2 - 7l e + 2A 2 ) 1 / 2 ' (33) 



kR-idk = 



V^(7i ~ 0A so 



£(2t; 2 fc 2 + 7l 2 + 7l £ + 2A 2 ) 1 /2' 



(34) 



and then assuming the appropriate limits of the inte- 
gration, one finds the final expressions for the spin Hall 
conductivity as presented below. 

When the chemical level is inside the gap, \fi\ < A so , 
the spin Hall conductivity is equal to 



<i = - V 



(35) 



When A so < | M | < v/ 7l 2 + A 2 D , 



2( 7l + v ^ALj Asq £ 
2 x /M 2 -A 2 + 7i H 4tt' 



while for > \J"f 2 + A 2 Q one finds 

2(/x 2 - A 2 J - 7l 2 4A so e 



■''!/ 



4(^ 2 -A 2 )- 7l 2 4tt 



(36) 



(37) 
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The spin Hall conductivity inside the gap is now twice 
as large as that in the case of a single-layer graphene. 
The general behavior of the conductivity with position 
of the Fermi level, shown in Fig. 3 by the solid dark (solid 
red) line, is qualitatively similar to that for a single-layer 
graphene, i.e., outside the gap the spin Hall conductivity 
tends to zero with increasing while it remains con- 
stant and quantized inside the gap. This behavior is rea- 
sonable as the spin Hall conductivity is due to spin-orbit 
coupling, which is the same for both atomic monolayers. 
Note, the formula derived correspond to one Dirac point, 
while the figures include contributions from both Dirac 
points. 

The corresponding spin Nernst conductivity is given 
by the following formulas: 



When A so < \n\ < yffi + A% 



6 [i 1 



-7?-2/i 2 - 



47l M 



37i 



4 7i Vm 2 -AL + V+7i 2 -4AL 



.(38) 



For > v/7i 2 + AL, 

y 8A*„ + 6AL7i 2 ±jtl 2 ( 8A 'q + + 8 ^ 

M 2 (4A2 D +7i 2 -4/i 2 ) 2 



(39) 



In turn, when fi is in the gap, < A so , the spin Nernst 
conductivity vanishes, a s x z y = 0. The signs — and + in the 
above formulas correspond to the case when the chemical 
potential is negative or positive, respectively. 

Variation of the spin Nernst conductivity with the 
chemical level is shown in Fig. 3 by the solid gray (solid 
blue) line. Similarly as in a single-layer graphene^, the 
spin Nernst conductivity vanishes for the Fermi level in- 
side the gap, when the system is in the insulating phase, 
and becomes nonzero for the Fermi level inside the va- 
lence or conduction bands, when the temperature gradi- 
ent generates a longitudinal charge current. Note, the 
spin Nernst conductivity becomes divergent as [i ap- 
proaches edges of the energy gap. 

Conductivity in the tight binding model can be ob- 
tained in a similar way, although the corresponding for- 
mulas are cumbersome and will not be presented here. 
Instead of this we present some numerical results, which 
in Fig. 3 are shown by the dotted lines. Note, the spin 
Hall as well as spin Nernst conductivities in the effective 
model coincide very well with the results obtained from 
the tight binding model. 



B. The case of V / 

Let us consider now the case of V ^ 0. The procedure 
presented above for the effective model with V = can 




0.05 0.1 0.15 0.2 0.25 0.3 
V/t [10" 3 ] 

FIG. 4. (color on-line) Zero-temperature spin Hall conductiv- 
ity as a function of the vertical voltage V and position of the 
Fermi level for the effective model described by the parame- 
ters given in Fig. 2. Contributions from both Dirac points are 
included. The vertical dashed line corresponds to V, where 
transition from the spin Hall insulator to conventional insu- 
lator takes place. 



be easily extended to a nonzero vertical bias, V ^ 0. 
The difference is that now the degeneracy of the bands 
is lifted and we have 8 different bands, n = 1 — 8, which 
have to be taken into account. Thus, instead of Eq.(27) 
wc have now 



He) 



X(e) 



rin=i[ £ -- E " + M + ^sign(e)] i 



(40) 



with adequate x{ £ )- Following the procedure described 
above for V = 0, one can derive the corresponding ana- 
lytical formula. These formula, however, will not be pre- 
sented here as they are rather cumbersome, so we present 
only numerical results. Moreover, since the results in the 
tight binding model coincide with those obtained with 
the effective Hamiltonian, as shown above, we restrict 
the analysis below to the effective Hamiltonian. 

In Fig. 4 we show the spin Hall conductivity as a func- 
tion of the Fermi level \i and vertical bias V . For V = 
we recover the quantized conductivity in the gap. As 
V increases, however, the range of quantized spin Hall 
conductivity shrinks and at a certain value of V (indi- 
cated by the dashed line in Fig. 4) there is a transition (at 
[i = 0) from a^ y = — 4(e/4-7r) to er£~ = 0. This behavior is 
explicitly shown in Fig. 5, where several cross-sections of 
Fig. 4 along constant values of \i are presented. The above 
transition is clearly evident for the curve corresponding 
to n = 0. 

The transition from topological insulating phase at 
small voltages to the normal insulating behavior at large 
voltages is also clearly visible in Fig. 6, which presents 
some cross-sections of Fig. 4 along constant values of V. 
This figure shows how the range of the quantized value of 
a^y changes with increasing V . As V increases starting 
from V = 0, width of the range where a^ z y is quantized 



FIG. 5. (color on-line) Zero-temperature spin Hall conductiv- 
ity as a function of vertical bias voltage V for indicated values 
of the Fermi energy. The curves correspond to crossections of 
Fig. 4 along constant ji. The other parameters as in Fig. 4. 



FIG. 7. (color on-line) Electronic spectra of the bilayer 
graphene in the vicinity of the Dirac point for different values 
of V . The curves show variation of the gap with increasing 
V, and closure of the gap at some critical value of V . The 
other parameters as in Fig. 4 
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FIG. 6. (color on-line) Zero-temperature spin Hall conductiv- 
ity as a function of the Fermi level for indicated values of the 
vertical bias voltage V. The curves correspond to crossections 
of Fig. 4 along constant V. The other parameters as in Fig. 4. 



FIG. 8. (color on-line) Low-temperature spin Nernst conduc- 
tivity as a function of the Fermi level for indicated values of 
V, calculated for the same parameters as in Fig. 6. 



shrinks, and at a certain critical value of V width of this 
range goes to zero. The spin Hall conductivity a s xy at 
li = changes then from o s x ~ y = —4(e/4n) at voltages 
smaller then the critical one to a s xy = at higher volt- 
ages. This clearly reveals a transition from the topolog- 
ical insulating phase to the normal insulating behavior 
(more information on the topological insulating phases 
in graphene can be found eg. in Ref. (32j). From Figs 
4 to 6 one could conclude that the gap diminishes with 



increasing V, becomes totally suppressed at the critical 
value of V , and then becomes open again at larger volt- 
ages. Indeed, this is the case as shown in Fig. 7, where 
the spectrum near the gap is plotted for several values of 
V. This figure clearly shows that the gap becomes closed 
at the critical value of V, and then is open again at larger 
values of V. The spin Hall conductivity in the gap above 
the critical voltage is however suppressed. 

The phase transition from the topological spin Hall 
insulating phase to the conventional insulator becomes 
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revealed in the spin Nernst conductivity, too. This is pre- 
sented in Fig. 8, where the low-temperature spin Nernst 
conductivity is shown as a function of the Fermi energy 
for the same values of V as in Fig. 6. The range of zero 
spin Nernst conductivity decreases with increasing V, 
goes to zero at the critical value of V, and then becomes 
nonzero again for larger values of V. 



IV. LOW-ENERGY EFFECTIVE 
HAMILTONIAN 

As we have already mentioned above, separation of 
the bands E% and £4 (or -E3.3' and Kt,4') is much larger 
than separation of the bands E\ and E2 (or E\ t y and 
£»2,2')- The latter determines the energy gap induced 
by the spin-orbit coupling (see Fig. 2). When only the 
electron states near the Fermi level are relevant and the 
Fermi level is in the gap or close to it, one may further 
reduce the effective Hamiltonian to include explicitly the 
bands E\^y and E2.2', and the other bands only via ef- 
fective parameters of the corresponding reduced effective 
model. The relevant reduced low-energy effective Hamil- 
tonian takes the form 2 ^: 



HI 



A so + V 




2m 








V 



2m 





-A, 



V 





h 2 k 2 _ 
2m 









A.,„ - V 



(41) 



where k± = k x ± ik y and m — ^i/2vf is the effective 
electron mass. In the following we consider some special 
cases. 



A. The case of V = 

The corresponding eigenvalues of the Hamiltonian (41) 
are then equal to 



Ei,2 — T 



A 



h 2 k 2 
2m 



1/2 



(42) 



see the inset in Fig. 9. Using the notation introduced in 
the preceding section we find 



de 



d 2 k 



where 



2 J 2nJ (2^) 2 n:;: i [e-i?„+M + ^sgn( £ )]4' 

(43) 



X(e) 



4A,„?i 4 fc 2 



in- 



h 2 k 2 

2m 



(44) 



Performing the integration over e and then over k, one 
arrives at the following analytical formulas for the spin 




FIG. 9. (color on-line) Spin Hall and spin Nernst conductiv- 
ities of a bilayer graphene for the low energy reduced Hamil- 
tonian with the parameters corresponding to those assumed 
in Fig. 2. Contributions from both K and K' points are in- 
cluded. The solid lines correspond to the results obtained 
from the Hamiltonian (10), while the dotted ones from the 
reduced Hamiltonian (41). The inset shows the energy spec- 
trum corresponding to Eq.(42). 



Hall conductivity: 



for > A so , and 



2A. S 



xy 



m 47r 



<*v = "IT 



(45) 



(46) 



when the chemical level is inside the gap, < A so . 

The corresponding low-temperature spin Nernst con- 
ductivity is then given as follows. For > A so one 
finds 



^6 



A s 



-T. 



(47) 



while for fi inside the gap a^ z y = 

The above results for both spin Hall and spin Nernst 
conductivities are shown by the dotted lines in Fig. 9, 
where they are compared with the corresponding results 
obtained from the effective Hamiltonian (10) (solid lines). 
There is a nice agreement between the results. We note 
that the divergency of the spin Nernst conductivity when 
the Fermi level approaches the band edges, observed in 
the tight-binding and effective Hamiltonian, is not repro- 
duced by the reduced Hamiltonian. 



B. The case of V / 

When V 7^ the eigenvalues of the Hamiltonian (41) 
take the form, 



E\a — T 



(V + A so ) 2 + 



h 2 k 2 

2m 



1/2 



(48) 
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Ey 



=F 



(^-A so ) 2 + 



h 2 k 2 

2m 



1/2 



(49) 



The nonzero V leads to splitting of the electron bands, 
as already discussed above. 

When > (V + A so ), the spin Hall conductivity is 
then given by 



s _ _ e 2A so 
axy ~ 4tt ]//| 

When | V - A ao | < \fx\ < (V + A so ), 

V - A. 



l/'l 



(50) 



(51) 



Finally, when [i is inside the gap, 

tr**=0 for l/>A so (52) 
a% = -1^ for F<A so . (53) 

The corresponding low-temperature spin Nernst con- 
ductivity is given by 



a xy — G K B 1 II 

y 6 n\n\ 
for > V + A so , and 



~k 2 B Tj±-{V- A so ) 



(54) 



(55) 



for | V — A. so | < < (V + A so ). In turn, the spin Nernst 
conductivity vanishes inside the energy gap. 

Behavior of spin Hall and spin Nernst conductivities 
with position of the Fermi level and bias voltage V al- 
most coincides with that presented in Figs 4 to 6 and 8, 
and therefore will not be present here. As in the k • p 
model, one observes the same transition between the spin 
Hall insulator and classical insulator as V increases. As 
already mentioned above, there is no divergency of the 
spin Nernst conductivity at the band edge for the low 
energy effective model considered here. 



SUMMARY 



We have calculated analytically as well as numerically 
the spin Hall and spin Nernst conductivities in a bilayer 
graphenc. To describe the relevant electronic spectrum 
we have assumed the tight binding model as well as some 
simplified effective Hamiltonians relevant for states close 
to the Dirac points. Both spin Hall and spin Nernst 
effects consist in transverse spin accumulation (spin cur- 
rent). However, as the spin Hall effect is due to exter- 
nal electric field, the Nernst effect is due to a tempera- 
ture gradient. Assuming intrinsic spin orbit interaction, 
we have found the intrinsic contributions to both effects. 
Generally, the spin Hall conductivity in the energy gap of 
a bilayer graphene is twice as large as that for a single- 
layer graphene. When external voltage is applied be- 
tween the two atomic sheets, we found a transition from 
the spin Hall insulator phase to the conventional insula- 
tor behavior. The energy gap at the transition point is 
closed. 

We have compared results obtained from the tight- 
binding model with those derived from the effective 
Hamiltonians. From this comparison we arrived at the 
conclusion, that both spin Hall and spin Nernst conduc- 
tivities are well described by the effective Hamiltonians, 
which allow derivation of some analytical results. 
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